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Abstract 


This  report  describes  a program  module  of  the  expert  system  EXPOSURE.  This  module,  written 
in  Common  Lisp,  is  for  calculating  the  maximum  electromagnetic  radiation  incident  upon  a 
building’s  wail  due  to  the  burning  of  a neighboring  building.  It  is  assumed  that  the  burning 
building  has  an  arbitrary  number  of  rectangular  openings  emitting  radiation.  The  exposed  wall 
can  be  considered  as  being  composed  of  a number  of  rectangular  regions.  These  regions  may  be 
openings  or  regions  of  different  materials.  This  program  can  determine  the  maximum  radiation 
for  each  of  these  regions  and  for  the  wall  less  these  regions. 


Keywords:  exposure,  radiation,  radiative  heat  transfer,  computer  programs. 


Forward 


In  this  report  we  are  trying  a somewhat  unusual  approach  to  documenting  a computer 
program.  This  entire  report  is  the  contents  of  the  computer  file  that  contains  the  program 
^’Radiation."  The  only  things  omitted  are  some  semicolons  that  indicate  that  the  wntten  texts 
of  this  report  are  comments  and  the  initial  header  strip: 

;;  Mode:  Joshua;  Package:  FIRE-PROTECTION;  Syntax:  Joshua;  Default-character- style: 
;;  (:FIX  :ROMAN  :LARGE)  -*-  Created  7/27/89  16:15:45  by  Smith  running  on  SMOKEY  at 
;;  NBS. 

The  only  addition  to  this  report  that  is  not  in  the  file  is  the  figure.  A computer  pro^am  is  a 
^namic  thing  and  should  be  documentated  during  development  as  well  as  when  it  is  finished. 
This  report  was  generated  by  printing  a file  which  is  a module  in  an  expert  system . It 
represents  a snap  shot  of  the  module  at  one  instance  of  time.  Since  the  computer  code  is  the 
focus  of  this  report  it  is  included  in  the  body  of  the  report.  The  document  is  designed  to  be 
read  in  its  entirety  although  the  audience  wanting  to  do  so  is  probably  small.  We  believe  this 
type  of  documentation  is  required  for  modules  of  expert  systems  which  are  going  to  be  used 
to  make  real  world  recommendations.  This  type  of  document  will  allow  others  to  review, 
revise,  and  maintain  this  code  and  its  physical  science  basis.  Those  who  do  not  understand 
LISP  code  can  skip  reading  the  code,  which  is  printed  in  bold  type  style,  enclosed  in 
parentheses,  and  often  starts  with  d«f . . . . These  readers  should  be  able  to  follow  the  logic 
of  the  program  but  they  may  not  be  able  to  determine  if  the  code  performs  as  described. 
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1.  Introduction 

This  report  documents  a module  of  the  expert  system  EXPOSURE.  It  describes  a program  for 
calculating  the  maximum  electromagnetic  radiation  incident  upon  a building’s  wall  due  to  the 
burning  of  a neighboring  building.  It  is  assumed  that  the  burning  (exposing)  building  has  an 
arbitrary  number  of  rectangular  openings  emitting  radiation.  Tlie  ex^sed  wall  can  be 
considered  as  being  composed  of  a number  of  rectangular  regions.  iTiese  regions  may  be 
openings  or  regions  of  different  materials.  We  wish  to  determine  the  maximum  radiation  for 
each  of  these  regions  and  for  the  wall  less  these  regions. 

After  rescuing  people  from  a burning  building,  the  next  highest  priority  item  for  fire  fighters 
is  to  prevent  tiie  fire  from  spreading  to  adjacent  buildings  [1].  Because  of  the  potential 
damage  that  can  be  caused  if  fire  propagates  from  buildmg  to  building  there  are  standards 
that  provide  guidance  on  how  to  reduce  the  likelihood  of  tire  spreading  [2]  [3]. 

A critical  part  of  the  analysis  of  the  spread  of  fire  between  buildings  is  the  calculation  of  the 
electromagnetic  radiation  falling  upon  the  exposed  surfaces  of  the  nonbuming  building.  The 
calculations  of  the  radiation  is  complex  and  does  not  lend  itself  to  simple  approximations  [5] 
[6]  [7]  [8].  Present  standards  are  based  on  work  done  before  the  advent  of  the  personal 
computer  [4]  [5].  To  make  the  estimating  of  the  radiation  level  in  the  standards  possible,  a 
number  of  simplifying  assumptions  had  to  be  made  in  the  existing  standards. 

NFPA  80A,  Recommended  Practice  for  Protection  of  Buildings  from  Exterior  Fire  Exposures 
[2],  makes  a number  of  assumptions.  Some  are  used  in  this  expen  system  others  are  not.  We 
do  make  the  following  assumptions  which  are  also  made  by  NFPA  80A: 

1.  Only  one  wall  of  the  burning  building  will  expose  one  wall  of  the  exposed  building. 

2.  All  radiating  openings  are  assumed  to  have  a rectangular  shape  [5]. 

We  do  not  make  the  following  assumptions  that  are  made  by  the  NFPA  80A. 

1.  The  exposed  and  exposing  walls  are  parallel  (for  an  exception  see  [3]). 

2.  A region  of  a wall  with  a number  of  radiators  can  be  treated  as  a single  radiator 
with  a size  of  the  rectangle  which  encloses  all  the  individual  openings  and  which 
radiates  at  a reduced  intensity.  The  reduced  intensity  equals  the  ratio  of  the  area  of 
the  openings  to  that  of  the  enclosing  rectangle. 

3.  If  the  exposed  building  has  a combustible  wall,  then  the  minimum  allowed  separation 
between  buildings  is  the  separation  such  that  the  thermal  radiation  falling  on  the 
exposed  building  just  equals  1.26  w./sq.  cm.  (the  piloted  ignition  level  for  oven  dried 
wood). 

4.  Only  the  maximum  radiation  level  on  the  exposed  wall  is  determined. 


Briefly,  this  program  does  the  following:  Based  upon  the  size  of  an  opening  and  its  minimum 
distance  from  the  exposed  wall,  the  program  divides  the  opening  into  subelements.  It  then 
finds  the  configuration  factor  due  to  one  opening  by  summing  the  contributions  from  all  the 
small  elements  of  the  area  in  the  opening.  It  then  sums  all  the  contributions  from  all  the 
openings.  This  gives  the  configuration  factor  at  any  point  on  the  exposed  wall  due  to  all  the 
openings.  The  program  then  can  determine  the  maximum  configuration  factor  in  any 
rectangular  region  on  the  exposed  wall  by  scanning  the  region.  Finally,  it  has  the  capability 
of  finding  the  maximum  configuration  on  the  area  of  the  exposed  wall  less  its  openings. 

The  program  described  in  this  report  was  written  in  Symbolics  Common  LISP.  (Note:  Certain 
commercial  products  are  identified  in  this  report  in  order  to  adequately  specify  the  computer 
program.  Such  identification  does  not  imply  recommendation  by  the  National  Institute  of 
Standards  and  Technology,  nor  does  it  imply  that  these  products  identified  are  necessarily 
the  best  available  for  the  purpose.) 
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2.  Source  Flux 

We  assume  all  surfaces  of  openings  in  the  exposing  building  radiate  according  to  the  equation 

E - zkT^  w/sq  cm.  ( i ) 

where  e is  the  emissivity  (e  < 1).  For  an  ideal  black  body  e = 1 and  where  X is  the  Stefan 
Boltzmann  constant  with  a value  of  5.6703e-12  (w./sq.  cm.)(K)'^-4.  We  will  assume  e = 1 for 
all  the  openings.  Since  a real  opening  will  always  radiate  less  than  a black  body  at  the  same 
temperature,  mis  assumption  results  m overestimating  the  radiation  and  therefore  errs  on  the 
side  of  safety. 

We  record  the  Stefan  Boltzmann  constant  in  the  program  as  follows: 

(cie£const:ant^  *St®fan-Boltzxaann-constant*  5.6703«-12  "(w/sq  cm)  (K)'^-4") 


The  LISP  function  biack-body-fiux  is  the  LISP  version  of  eq.  1.  It  takes  as  its  input  or 
its  parameter  the  temperature  of  the  burning  compartment.  It  returns  the  radiation  flux  at  the 
surface  of  the  openings  of  the  building  in  w./sq.  cm. 

(daftm  black -body -flux  (bamporabura)  Temperature  ifi  Kelvin 

*Sbafan-Bolb2maxm-consbanb*  (expb  bamparabura  4} ) ) 


3.  Configuration  Factor 

If  we  assume  Lambert’s  cosine  law  holds,  the  incident  radiant  flux  falling  on  a small  element 
of  surface  is 

Rf  =»  E j (cos  CLl  cos  a2)  / ( n r''2)  dAl  (2) 

where  we  integrate  over  the  surface  ai,  and  e is  the  source  flux  determined  above.  The 
integral  is  called  the  configuration  factor.  From  Figure  1.  we  see  that  r is  the  distance 
between  the  point  of  integration  at  dAi  to  the  point  at  which  we  are  evaluating  the  field,  ai 
is  the  angle  between  the  normal  to  dAi  (the  exposing  surface)  and  r,  and  a2  is  the  angle 
between  the  normal  to  the  exposed  surface  at  tne  point  we  are  evaluating  the  field  and  r.  If 
(xa,  ya,  za)  is  the  source  point  on  dAi  and  (xO,  yO,  zO)  is  the  field  point  on  the  exposed 
surface,  then  we  have 

r = ( (xO  - xa)''2  + (yO  - ya)  ^2  + (zO  - za)^2)^.5 

If  we  define  y and  pi  such  that 


cosy  “ 

(xO  - 

■xa)  /r 

slny  = 

(yO  - 

■ ya)/r 

cosPi 

= (x2 

- Xl)/Ll 

sinPi 

= (y2 

- yl)/Ll 

where  pi  is  the  slope  of  walll  (the  exposing  wall)  and  the  length  of  the  wall  is  given  by 
Li  = ((x2  - xl)''2  + (y2  - yl)'^2)''.5 
then  we  have 


cos  ai  = cosy  sln\il  - siny  cosPi  . 
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Figure  1.0  Variables  of  the  Configuration  Integral 
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If  we  define  P2  such  that 

cosp2  = (x4  - x3) /L2 
slnP^  - (y4  - y3)/L2 

where  P2  is  the  slope  of  wall2  and  the  length  of  exposed  wall  (wall2)  is  given  by 

L2  - ( (x4  - x3)''2  + (y4  - y3)''2)''.5 

then  we  have 

cos  0.2  - cosy  s±n^2  - slay  cos^2 . 

All  rectangles  are  defined  by  giving  the  coordinates  of  the  lower  left  hand  comer  and  the 
upper  right  hand  comer.  This  is  sufficient  to  define  a rectangle  because  all  rectangles  of 
interest  are  either  vertical  or  horizontal.  Therefore,  an  opening  or  a wall  can  be  defined  by 
the  four  points:  (xi,  yi,  zi),  (x2,  y2,  zi),  (x2,  y2,  z2;,and  (xi,  yi,  z2)  or 
by  using  the  two  point  convention,  (xi^  yi,  zi) 
and  ix2,  y2,  z2). 


3. 1 Program  Input  and  Basic  Functions 

The  input  for  this  program  is  the  lower  left  and  upper  right  comers  of  the  exposing  wall  and 
its  openings  and  the  exposed  wall  and  its  openings.  Important  concepts  for  this  program  are 
various  types  of  points  and  rectangles.  In  this  section  we  give  their  definitions  and  the  basic 
functions  for  operating  on  them. 

One  basic  entity  is  the  point  which  we  define  as: 


(dofserucb  (polne  :al'tarane) 
at  y z cf) 

Sometimes  we  have  an  interest  in  the  configuration  factor  for  a point  so  it  has  been  included. 
Another  basic  entity  is  the  vertical  rectangle  which  we  define  as: 

(defstaruct:  (vaarblcal-racbangla) 
lowar-lafe  uppar-rlght:) 

Then  walls,  openings,  and  elements  of  integration  are  just  instances  of  vertical  rectangles. 

(daf  struct,  (wall  ( : Includa  vaartlcal-ractangla)  ) ) 

(daf struct  (opanlng  (: Includa  vartlcal-ractangla) ) ) 

(daf struct  (alamant  : altarant  (: Includa  vartlcal-ractangla))) 


We  will  have  use  for  two  simple  functions.  One  is  the  squara-root  and  the  other  is  the 
squara.  We  define  these  functions  as  follows: 

(defsubst  squara-root  (n\unbar) 

(sqrt  numbar) ) 

(dafsubst  squara  (numbar) 

(axpt  numbar  2) ) 

We  used  dafsubst  SO  they  get  expanded  at  compile  time. 
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We  will  have  use  for  a function  which  rounds  off  some  of  our  results  to  a controllable 
number  of  significant  figures.  We  call  that  function  rotmd-of  f . 


(dafuxi  ro\md-o££  (numb«x’  £op^i.onaJ.  (si.gxxi.fi.cane-flguz’as  0) ) 

(lab  ( (d  (azpe  10  slgiLiflcane-flguras) ) ) 

(/  (round  (*  numbar  d) ) d) ) ) 

We  will  have  need  to  determine  the  separation  between  two  points. 

(dafun  polnb-saparablon  (poinb-l  polja'b-2)  ,'3  dimensions 
(squara=roob 

(+  (squara  (-  (polnt-x  polnt:-2)  (polnb-x  poinb-l) ) ) 

(squara  (-  (point-y  point-2)  (point-y  point-1))) 

(aquara  (-  (point-z  point-2)  (point-z  point-1)))))) 

We  use  point-saparation  to  determine  the  length  of  the  diagonal  of  vertical  rectangles. 

(dafun  vartical-ractangla-diagonal  (vartical-ractangla) 

(point-saparation  (vartical-ractangla-lowar-laft  vartical-ractangla) 

(vartical-ractangla-uppar-right  vartical-ractangla) ) ) 

We  will  also  need  the  separation  of  the  projection  of  two  points  on  the  xy-plane. 

(dafun  xy-saparation  (point-1  point-2)  ,'SepCirCltion  in  Xy  plane 
( squara-root 

(+  (squara  (-  (point-x  point-2)  (point-x  point-1))) 

(squara  (-  (point-y  point-2)  (point-y  point-1)))))) 

We  use  xy-saparation  to  determine  the  width  of  vertical  rectangle. 

(dafun  vartical-ractangla-width  (vartical-ractangla) 

(xy-saparation  (vartical-ractangla-lowar-laft  vartical-ractangla) 

(vartical-ractangla-uppar-right  vartical-ractemgla) ) ) 

Instead  of  calling  the  wall  dimension  in  the  xy  plane  its  width,  we  called  it  length. 

(dafun  wall-langth  (wall) 

(vartical-ractangla-width  wall) ) ;nOte  Switchform  length  tO  width 

But  for  an  opening  we  switch  back  to  width. 

(dafun  opening-width  (opening) 

(vartical-ractangla-width  opening) ) 

We  will  need  the  height  of  our  rectangles  so  we  define  the  following: 

(dafun  vartical-ractangla-haight  (vartical-ractangla) 

(-  (point-z  (vartical-ractangla-uppar-right  vartical-ractangla)) 

(point-z  (vartical-ractangla-lowar-laft  vartical-ractangla)))) 

(dafun  wall-height  (wall) 

(vartical-ractangla-haight  wall) ) 

(dafun  opening-height  (opening) 

(vartical-ractangla-haight  opening) ) 

Likewise,  we  will  need  the  area  of  openings  and  elements  of  integration. 

(defun  vartical-rectetngle-srea  (vertical-rectangle) 

(*  (vertical-rectangle-haight  vertical-rectangle) 
(vertical-rectangle-width  vertical -rectangle) ) ) 
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(dafun  op«n±ng-a]r«a  (op«iilAg) 

(vairt:i.oaX-r«ct:angla-araa  opanijig)  ) 

(dafun  alamant-araa  (alamant) 

(vart±C2Ll~rac^amgla-araa  aXamanfe) ) 

The  following  functions  are  used  to  calculate  the  coordinates  of  the  centers  of  various 
rectangles. 

(dafun  var^±caX-ract:angla-cant,ar  (vairtlcal-rac^angla) 

(la^  ( (11  (vaz-klcal-racbangla-lowar-laft.  vasrtleal-rac^angla) ) 

(ur  (vairblcal-ractangla-uppar-rlght  vairklcal-rac^angla)  ) ) 
(maka-polnfc  :ac  (/  (+  (polnt~x  ur)  (poiat-x  11))  2) 

•y  (/  (+  (polnt-y  ur)  (point-y  11))  2) 

IX  (/  (+  (point-*  ur)  (point-*  11))  2)))) 


(dafun  opaning-cantar  (opaning) 

(vartical-ractangla-cantar  opaning) ) 

(dafun  alamant- canter  (alamant) 

(vartical-ractangla-cantar  alamant) ) 

The  slope  of  the  various  rectangles  in  the  xy-plane  can  be  determined  by  these  functions. 

(dafun  slopa-of-vartical-ractangla  (vertical -rectangle) 

(let*  ( (11  (vartical-ractangla-lowar-laft  vertical-rectangle) ) 

(ur  (vartical-ractangla-uppar-right  vertical -rectangle) ) 

(y2  (point-y  ur) ) 

(yl  (point-y  11) ) 

(x2  (point-x  ur) ) 

(xl  (point-x  11) ) ) 

(atan  (-  y2  yl)  (-  x2  xl)  ) ) ) 

(dafim  slopa-of-wall  (wall) 

(alopa-of-vartical-ractangla  wall)) 

(defiin  alopa-of -opaning  (opaning) 

(slopa-o£-vartical-ractangla  opaning) ) 

(daf\m  slope -o£-alamant  (alamant) 

(slopa-o£-vartical-ractangla  alamant) ) 

This  function  determines  the  cosines  of  the  alpha  angles. 

(dafun  cos-alphas  (alamant  field-point  field-wall) 

(let*  ( (al-cantar  (alamant -canter  alamant) ) 

(r  (point- separation  al-cantar  field-point) ) 

(cosy  (/  (-  (point-x  field-point) 

(point-x  al-cantar) ) r) ) 

(siny  (/  (-  (point-y  field-point) 

(point-y  al-cantar) ) r) ) 

(Pl  (slopa-of-alamant  alamant) ) 

(P2  (slopa-of-wall  field-wall) ) 

(cos-cxl  (-  (*  cosy  (sin  pi)  ) 

(*  siny  (cos  pi)))) 

(cos-a2  (-  (*  cosy  (sin  P2)  ) 

(*  siny  (cos  p2))))) 

(values  (abs  cos-al)  (abs  cos-a2) ) ) ) 
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3.2  Configuration  Factor  of  an  Element 

Whenever  we  have  the  value  of  the  diagonal,  d,  of  an  element  of  an  opening  much  less  than 
the  value  of  the  shortest  value  of  r over  the  opening,  i.e.,  d « r-min  where 

d - ( (x2  - xl)^2  + (y2  - yl)''2  + (z2  - zl)'^2)'',5 

then  ai , (x2,  and  r do  not  vary  appreciably  over  the  element  of  area  ai. 

Therefore,  we  can  approximate  the  configuration  integral  for  one  element  by  the  following 
expression: 


cf  - (cos  0.1  cos  02)  Al/  (pi  r''2) 


We  will  define  p±  to  be  3.141593.  This  avoids  double  precision  calculations  and  thereby 
speeds  up  the  computations. 

(dafvax  *singl«-p-pl*  3.141593) 

We  now  define  the  LISP  function  configuration-factor-for-ona-elomant  with 
parameters  of  aiamant,  f i«id-point,  and  fiaid-waii.  It  returns  the  numerical  value 
for  the  configuration  factor  for  this  infinitesimal  element  at  the 

fi.ald-poi.ne. 

(dafun  configuration- factor- for- ona-alamant  (alamant  fiald-point  field-wall) 
(multipla-vetlua-bind  (coaal  coao2) 

(cos-alphaa  alamant  fiald-point  fiald-wall) 

(/  (*  cosal  cosa2 

(alamant-araa  alamant) ) 

(*  (square  (point- separation  (alamant-cantar  alamant)  fiald-point) ) 
*singla-p-pi*) ) ) ) 

We  have  found  that  if  we  have  d < .1  * r-min,  we  do  not  need  to  integrate  over  the  area 
with  this  value  of  d.  We  get  a good  approximation  for  our  integral  by  using 
configuration-factor-f or-ona-aiamant.  Therefore  we  define  our  constant 
*limit-alamant-diagonal-to-r-min-ratio*  aS  folloWS 

(dafvar  *limit-elamant-diagonal-to-r-min-ratio*  .1) 


If  our  area  of  interest  does  not  satisfy  this  condition,  we  must  subdivide  it  into  subelements 
until  it  does. 

To  implement  this  approximation  we  need  to  be  able  to  1.)  calculate  the  shortest  distance 
between  a point  and  a planar  region  (r-min)  and  2.)  determine  the  number  of  subelements 
so  we  can  use  the  above  approximation  to  calculate  the  configuration  factor. 


3.3  Shortest  Distance  Point  to  Planar  Region 

We  wish  to  determine  the  shortest  distance  between  a point  (xo  yo  zO)  and  a planar 
rectangular  region  defined  by  its  coordinates  of  opposite  comers  (xi  yi  zi)  and  (x2  y2 
z2).  We  assume  the  region  lies  in  a plane  that  is  perpendicular  to  the  xy  plane,  i.e.  the 
vertical  plane. 


The  inputs  for  the  function  short«st-dlstanc«-point-to-plama-region  are  the  point 
and  varticai-ractangia.  It  retums  the  shortest  distance  and  the 
coordinates  of  the  nearest  point  in  the  rectangle. 


(defun  short es t -distance -point -to-plane- region 
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;;for  vertical  rectangle  in  rectangular  world 

(po±n't  vaftlcal-rac^angla) 

(i«t*  ( (xi  0)  :xi  etc.  will  be  the  nearest  point  in  the  region. 

(yi  0)  ;Value  ofO  for  xi,  yi,  & zi  are  to  avoid  nil. 

(zi  0) 

(zll  (point-z 

(v«irt:±cal-racbangltt-low«r-lo£fc  vartlcal-roctiangla) ) ) 

(zu  (poiat-z 

('^«zrki.cal-rac£angla-upp«r-r±gh£  vaxtlcal-rac^anglo) ) ) 

(yl  (point-y 

(var^lcal-roctangla-lowar-laft:  varlii.cal-rac'bangla) ) ) 

(y2  (point-y 

(varticai-ractangl«-upp«r-right  vartlcal-ractangla) ) ) 

(yll  (min  yl  y2) ) 

(yu  (max  yl  y2) ) 

(xl  (point~x 

(vaxticai-ractangla-lowar-laft  vartical-ractangla) ) ) 

(x2  (point-x 

('7’articai-ractangla-uppar-right  vartical-ractangla)  ) ) 

(xll  (min  xl  x2)) 

(xn  (max  xl  x2 ) ) 

(yO  (point-y  point) ) 

(xO  (point-x  point) ) ) 

;;  This  condition  on  z works  only  for  vertical  walls 

(aatq  zi  (min  zu  (max  zll  (point-z  point) ) ) ) 

;;  The  X and  y coordinates  of  a point  in  the  planar  region  are  not  independent. 

;;  They  are  related  by  the  equation  y = mx  b where  m is  the  slope  given 
;;  by  (y2  - yl  )l(x2  - xl ).  We  need  to  deal  with  two  special  cases  whim  are 
;;  when  m = 0 and  llm  = 0. 

(cond  ((-  x2  xl)  .‘Vertical  slope=  oi  because  xi  = xl . 

(aatq  xi  xl)  .’All  x’ s have  some  value  on  plane  of  region. 

(aatq  yi  (min  yu  (max  yll  yO)))) 

( (=  y2  yl)  ;m  = 0 SO  Urn  = <» 

(aatq  yi  yl) 

(aatq  xi  (min  xu  (max  xll  xO) ) ) ) 

;;Now  we  deal  with  the  general  case. 

(t  (lat*  ((m  (/  (-  y2  yl)  (-  x2  xl) ) ) .’Slope  of  Opening  in xy  plane. 

;;  Substituting  yl  for  y and  xl  for  x and  solving  for  b we  have 

(b  (-  yl  (*  m xl))) 

;;The  equation  for  the  line  which  is 
.’.perpendicular  to  above  line  is  y = (llm)x+  bl. 

.’.’Since  the  point  (xO,  yO)  lies  on  this  line 
;;  we  have 

(bl  (+  yO  (/  xO  m))) 

;;The  intersection  of  these  two  lines  is  given  by 

(xii  (/  (-  bl  b)  (+  m (/  1 m)))) 

(yii  (/  (+  b (*  m m bl) ) (+  (*  m m)  1)))) 

(cond  ( (and  (or  (<  xll  xii)  (=  xll  xii) ) 

;;The  intersection  is  in  the  region 

(or  (<  xii  xu)  (=  xu  xii) ) ) 

(aatq  xi  xii) (aatq  yi  yii) ) 

( (<  xii  xll) (aatq  xi  xll) (aatq  yi  (+  (*  m xi)  b) ) ) 

( (>  xii  xu) (aatq  xi  xu)  (aatq  yi  (+  (*  m xi)  b) ) ) ) ) ) ) 

(liat 

(aquara-root  (+  (aquare  (-  xO  xi) ) (aquare  (-  yO  yi) ) 

(aquara  (-  (point-z  point)  zi) ) ) ) 

(maka-point  :,x  xi  :y  yi  :z  zi)  ) ) ) ;nearest  point 
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3.4  Number  of  Subelements  in  Opening 

We  have  an  opening  with  sides  longer,  l,  and  shorter,  s,  we  want  the  diagonal  to  be  less 
than  *lljniit-element-dlagonal-to-r-min-ratlo*  times  the  minimum  V^Ue  of  the 
distance  from  the  field  point  and  any  point  in  the  opening.  If 

( (L''2  + s''2)''.5)  < <'*lijait-eleiaant-eLlagonal-to-r-iiiln-ratio*; 

then  we  are  done.  However,  if  this  is  not  the  case,  we  need  to  subdivide  the  opening  into 
approximate  square  elements.  Therefore,  we  let  n - l/s,  then  if  we  divide  l by  n we  get 
square  elements.  If  we  divide  i,  by  w - (round  n;  we  get  nearly  square  elements.  If 

((L/n)''2  + s''2)''.5  < ^'*lijalt-«l®mont-di.agonal-to-r-m±n-ratio* ^ *r-raln, 

we  are  done.  Otherwise,  we  must  continue  making  our  elements  smaller  by  reducing  the  size 
of  both  dimensions  by  dividing  them  by  p. 

Then  we  want  an  integer  p such  that 

( ( (L/N) ''2  + s^2)^.5)/P  < 

^'*liinit-al«m«nt-d±agonal-to-r~min-ratio*;  *r-raln  < ( ( (L/N) ''2  + 
s''2)  ''.5)/(P  -1)  . 

This  can  be  achieved  as  long  as  r-min  is  not  zero.  If  r-min  is  small,  p may  become  so  large 
as  to  make  this  approach  undesirable. 

We  make 

p = ( ( (L/N)  ^2  + 5^  / ^*lljait-«l«n«nt-d±agonal-to-r-iain-ratio*^  . 

It  follows  that  p - (calling  p) . We  rote  that  (calling  p)  is  a LISP  expression  which 
returns  an  integer  which  equab  p if  p is  an  integer  otherwise  it  returns  the  first  integer 
larger  than  p. 

We  recall  that  LISP  can’t  tell  the  different  between  n and  n. 

The  function  subaiomants-parajnatars  takes  the  width  and  height  of  the  opening  and 
r-min  as  input  and  returns  a list  which  has  the  form: 

(divisor-of-width  divisor-of-height). 


(dafun  subalomanes-paraxnaears  (opanlng-wldeh  opanlng-halghe  r-mln) 

(cond  ( (=  r-min  0) 

(arror 

"I  can  noe  con^uea  tha  solution  for  fwo  touching  walls, 
i.a.,  r-min  = 0.")) 

(t  (lat*  ( (longar  (max  opaning -width  opaning-haight  ) ) 

(shortar  (min  opaning-width  opaning-haight) ) 

(N  (round  (/  longar  shortar) ) ) 

(P  (calling 

(/  (squara-root  (+  (squara  (/  longar  N) ) 

(squara  shortar) ) ) 

(*  *limit-alament-diagonal-to-r-min-ratio*  r-min) ) ) ) ) 
(cond  ( (>  opaning-width  opaning-haight) 

(valuas  (*  p N)  p) ) ;dmsor-of -width  divisor-of-height 

(t  (valuas  P (*  P N)))))))) 
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3.5  Configuration  Factor  for  One  Opening 


To  obtain  the  configuration  factor  for  one  opening,  we  sum  the  configuration  factor  for  all 
the  subelements  in  that  opening.  The  function  configuration- factor- for  - on® - opaning 
take  as  input  the  fi®id-point,  opaning,  and  axpoaad-waii.  It  returns  the  configuration 
factor  of  the  opening  at  the  f iaid-point. 


(dafun  configuration-f  actor-for-ona-opaning  ^On  exposcd  wall  at  field  point 
(fiald-point  opaning  axposad-wail) 

(lat  ( (width-of-opaning  (opaning -width  opaning) ) 

(haight-of-opaning  (opaning-haight  opaning) ) 

(Pl  (slop® -of -opaning  opaning) ) 

;;  rm  is  the  shortest  distance  from  the  opening  to  the  field  point. 

(rm  (car 

(shortaat-diatanca-point-to-piana-ragion 
fiald-point  opaning) ) ) ) 

(multipla-valua-bind  (diviaor-of-width  diviaor-of-haight) 

(aubalamanta-paramatara  width-of -opaning  haight-of-opaning  rm) 

(lat*  ( (alamant-langth  (/  width-of-opaning  diviaor-of-width) ) 

( alamant - haigh t 

(/  haight-of-opaning  diviaor-of-haight) ) 

;;height‘Of -opening  = (-  ^ zl) 

(xl  (point-x  (opaning-lowar-laft  opaning) ) ) 

(yl  (point-y  (opaning-lowar-laft  opaning) ) ) 

(zl  (point-z  (opaning-lowar-laft  opaning))) 

( lowar- la  ft -alamant -point 

(maka-point  ;x  xl  ;y  yl  :z  zl)) 

( uppar- right -alamant-point 

(maka-point  :x  (+  xl  (*  alamant-langth  (cos  |3l)  ) ) 

:y  (+  yl  (*  alamant-langth  (sin  pi) ) ) 

:z  (+  zl  (*  alamant -haight) )) ) 

(acanning-alamant  (maka-alamant  :lowar-laft  lowar-laft-alemant-point 

: uppar-right  uppar-right-alamant-poin 


t))) 


(loop  for  i from  0 to  (-  diviaor-of-haight  1) 

;,ioop  over  elements  in  opening 

with  aum  = 0 
do 

(altar-point  lowar-laft-alamant-point  z (+  zl  (*  i alamant -haight) ) ) 
(altar-point  uppar-right-alamant-point  z (+  zl  (*  (+  i 1)  alamant-haight ) ) ) 

(loop  for  j from  0 to  (-  diviaor-of-width  1) 
do 

(altar-point  lowar-laft-alamant-point  x (+  xl  (*  j alamant-langth  (cos  pi)  ) 


y (+  yl  (*  j alamant-langth  (ain  pi) ) ) ) 

(altar-point  uppar-right-alamant-point  x (+  xl  (*  (+  j 1)  alamant-langth  (c 

os  Pl))) 

y (+  yl  (*  (+  j 1)  alamant-langth  (sin  Pl) ) ) ) 

(altar-alamant  acanning-alamant  lowar-laft  lowar-laft-alamant-point 
uppar - right  uppar-right - alamant -point ) 

(satq  aum  (+  sum 

(configuration-factor-for-ona -alamant 

acanning-alamant  fiald-point  axposad-wall) ) ) ) 
finally  (ratum  sum)  ) ) ) ) ) 


Comparing  the  results  of  configuration-factor-for-ona-opaning  with  analytical 
expressions  for  three  special  cases  showed  agreement  to  4 or  5 decimal  places  (see  Appendix 
a).  Note  that  this  is  an  approximate  integration  which  depends  upon  such  things  as  the  value 
of  *limit-alamant-diagonal-to-r-min-ratio*.  The  exaCtneSS  of  this  approximation 
can  be  improved  by  the  proper  adjustment  of  this  parameter. 
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3.6  Configuration  Factor  for  N Openings 

If  we  have  an  arbitrary  number  of  openings,  we  sum  the  configuration  factors  for  all  the 
individual  openings.  The  function  configuration-factor-for-n-op^nings  has  input  of 
the  fiaid-point,  -op«n±ngs , and  axposad-waii  It  retums  the  Configuration 

factor  due  to  all  the  openings  at  the  field  point 


(dafun  conflgura'bloii-faobor-for-n-opanl.ngs 

(fJ.ald-poi.nt:  lls'b-of-opani.nga  axpoaad-wall) 

(loop  for  opanlng  Jja  llae-of-opanlnga 
simmlng 

(eonfiguration-f  actor-for-ona-opanljig 
flald-poljit  opanlng  azposad-wall) ) ) 

At  this  point  we  can  calculate  the  field  at  any  point  on  the  exposed  wall. 


4.  Determining  Location  and  Strength  of  Maximum  of  an 
Opening  on  an  Exposed  Wall 

The  next  st^  is  to  determine  the  locations  and  strengths  of  the  maximum  configuration 
factor  (or  field)  on  an  exposed  wall  due  to  the  radiation  from  one  opening  in  the  exposing 
wall. 

If  we  have  a group  of  openings  for  which  the  separation  between  openings  is  much  smaller 
than  the  shortest  distance  to  the  exposed  wall,  then  such  a block  of  openings  can  be  treated 
as  a single  opening  whose  area  is  the  sum  of  the  individual  opening’s  areas. 

It  can  be  shown  that  the  maximum  for  an  element  of  an  opening  will  lie  on  the  exposed  wall 
at  the  point  closest  to  the  center  of  the  element. 

To  find  the  the  maximum  of  an  opening  we  will  start  by  guessing  its  location.  Then  we  will 
search  for  a point  that  has  a configuration  factor  that  is  close  to  the  maximum  configuration 
factor.  Close  means  within  2%  or  less.  So  we  need  to  pick  a point  for  our  initial  guess  and 
the  size  of  steps  we  wiU  take  while  searching  for  our  maximum. 

Because  of  symmetry  in  the  z direction,  the  maximum  for  an  opening  will  always  have  a z 
value  equal  to  the  midpoint  height  of  the  opening  if  it  falls  in  the  exposed  rectangle  of 
interest.  In  addition,  the  maximum  of  an  opening  frequently  lies  at  the  point  which  is  nearest 
the  center  of  the  opening.  This  tends  to  be  a better  guess  the  smaller  the  alpha  angle. 
However,  it  is  still  reasonable  even  out  to  45  degrees.  Our  first  guess  for  the  location  for 
the  maximum  is  at  the  point  nearest  the  center  of  the  opening.  Then  we  try  to  walk  to  a 
greater  value  for  the  configuration  factor.  We  move  in  a plane  parallel  to  the  xy  plane  a 
distance  given  by  the  step  size  and  see  what  the  configuration  factor  is  at  that  point. 

The  function  st«p-siz«  takes  as  input  the  opening  and  the  «xpos«d-waii.  To  estimate 
the  step- sir®  we  assume  that  the  relative  change  in  the  configuration  factor  goes  as 
(atap- a Lz9lseparation)'^2  which  is  characteristic  of  an  inverse  square  relationship  where 
the  change  is  approximately  perpenticular  to  r.  We  believe  this  is  a conservative  assumption. 
Significant  improvement  in 

the  code’s  speed  may  be  achieved  by  an  improvement  here. 


(daf\in  step-size  (opening  exposed-wall) 

(let*  ( (f  .02)  :max  fractional  change  in  cf  2% 

(source-point  (opening- center  opening) ) 

(r  (first  (shortest-dlstance-polnt-to-plane-reglon 
source-point  exposed-wall) ) ) ) 

(*  r (square-root  f ) ) ) ) 


The  function  Maxlmum-of-an-Openlng  takes  aS  input  the  opening  and  exposed-wall 
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and  returns  the  f i«id-point  (the  location  of  the  maximum  for  the  opening)  with  the  cf 
slot  having  a value  of  cfmo  (the  maximum  configuration  factor  for  tne  opening). 

(d«f\in  Maxi  imim-of -an-op«n±ng  ;returns  the  mox-field point 

(opttnljig  axposod-wall) 

(l®t*  ( (can'tox-o£-op«iLlng  (opanlng-cant^ttx  op«ni.ng) ) 

(shoirei  (8hozre«st:-dl,s^anc«-polAt:-to-pXan«-ragi.on 
c«n^ ax - o £ - opanXng 
axposad-wall) ) 

(£Xald-poi.ni:  (sacond  ahorfcl)  ) ^nearest point 

(P2  (slopa-o£-wall  axposad-wall) ) 

(stap  (s£ap-s±za  opani.ng  axposad-wall) ) 

(c£mo 

( configuration- f actor- for- ona - opanlng 
flald-polnt  opanlng  axposad-wall) ) 

(naw-flald-polnt+  (maka-naw-flald-polnt-posltlva-stap 

axposad-wall  flald-polnt  stap  |32)  ) 
(naw-flald-polnt-  (maka-naw-flald-polnt-nagatlva-stap 

axposad-wall  flald-polnt  stap  ^2) ) 

(cf+  (conflguratlon-factor-for-ona-opanlng 

naw-£lald-polnt+  opanlng  axposad-wall) ) 

(cf-  (conflguratlon-factor-for-ona-opanlng 

naw-flald-polnt-  opanlng  axposad-watll)  ) ) 

(cond  ( (>  cf+  (+  cfmo  .00001)) 

(saarch-for-max-plus 

axposad-wall  naw-flald-polnt+  stap  (32  cfmo  opanlng)  ) 

( (>  cf-  (-  cfmo  .00001))  ;This  and  cf+  make  the  following 

; ;approximate  equal  to, 

(saarch-for-max-mlnus 

axposad-wall  naw-flald-polnt-  stap  (32  cfmo  opanlng)  ) 

(t  (altar-point  flald-polnt  ;We  odd  info  to  the  field  point. 

cf  cfmo ) 

flald-polnt) ) ) ) ;No  larger  cf,  return  altered  field  point. 

We  first  try  to  move  in  a positive  x or  y direction.  If  we  don’t  find  a maximum,  we  look  in 
the  negative  direction.  We  note  that  for  one  opening  the  configuration  factor  is  monotonic  to 
its  maximum.  We  take  a step  by  making  a new  field  point  either  in  the  positive  or  negative 
direction.  We  first  write  the  function  for  a positive  step. 

(dafun  maka-naw-flald-polnt-posltlva-stap 
(axposad-wall  flald-polnt  stap  (32) 

(let  ( (x3  (polnt-x  (wall-lowar-laft  axposad-wall) ) ) 

(z4  (polnt-x  (wall-uppar- right  axposad-wall) ) ) 

(y3  (polnt-y  (wall-lowar-laft  axposad-wall) ) ) 

(y4  (polnt-y  (wall-uppar- right  axposad-wall) ) ) 

(xm  (polnt-x  flald-polnt) ) 

(ym  (polnt-y  flald-polnt) ) 

(zm  (polnt-z  flald-polnt) ) ) 

(cond  ( (<  (wall-langth  axposad-wall)  stap) 

flald-polnt)  ;No  new  field  points. 

( (not  (=  x4  x3 ) ) 

;;Dont  want  to  divide  by  zero  in  slope  calculation. 

(lat*  ( (slopa  (slopa-of-wall  axposad-wall) ) 

(b  (-  y3  (*  slopa  x3) ) ) 

;;  We  will  not  change  zm  but  only  xm  and  ym. 

(xpc  (+  xm  (*  (cos  (32)  (*  1 stap)))) 

(ypc  (+  (*  slopa  xpc)  b) ) ) 

;;  If  we  can  take  a step  in  the  positive  direction 
;;  and  stay  on  the  wall,  we  do  so. 

(cond  ((and  (or  (>  xpc  (min  x3  x4) ) (=  xpc  (min  x3  x4))) 

(or  (<  xpc  (max  z3  x4) ) (=  xpc  (max  x3  x4) ) ) ) 

(maka -point  :x  xpc  :y  ypc  ;z  zm) ) 
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(t  flald-poiai-b) ) ) ) 

(t  (i«t  ((xpcx4)  :For  vertical  lines,  all  x' s are  the  same. 

(ypc  (+  ym  (*  1 stap) ) ) ) 

;;  Try  to  walk  in  the  positive  y direction  on  the  wall 

(cond  ( (and  (or  (>  ypc  (mi.n  y3  y4) ) (=  ypc  (min  y3  y4) ) ) 

(or  (<  ypc  (max  y3  y4) ) (=  ypc  (max  y3  y4) ) ) ) 
(maka-polnt  :x  xpc  :y  ypc  :z  zm) ) 

(t  fiald-point) ) ) ) ) ) ) 

A negative  step  is  just  a positive  step  with  a negative  step  size. 

(dafun  maka-naw-flald-polnt-nagatlva-stap 
(axpo  sad- wall  flald-polnt.  stap  ^2) 
(maka-naw-fiald-polnt-positiva-stap 
axposad-wall  fiald-point  (*  -1  stap)  [32)  ) 


The  following  two  functions  are  used  above  to  search  either  in  the  positive  or  negative 
directions  for  the  maximum. 

(dafun  saarch-for-max-plus 

(axposad-wall  fiald-point  stap-slza  P2  cfmo  opanlng) 

(lat*  ( (naw-fiald-point+  (maka-naw-f iald-polnt-positiva-stap 

axposad-wall  fiald-point  stap-slza  (32)) 
(cf+  (configuration-factor-for-ona-opaning 

naw- fiald-point -t-  op>anlng  axposad-wall)  ) ) 

(cond  ( (>  cf+  (+  cfmo  .00001)) 

(saarch-for-max-plus 
axposad-wall  naw-flald-polnt+ 
stap-slza  {32  cf+  opanlng)  ) 

(t  (altar-point  fiald-point  ;We  add  info  to  the  field  point. 

cf  cfmo) 
fiald-point) ) ) ) 

(dafun  saarch-for-max-mlnus 

(axposad-wall  fiald-point  stap-slza  (32  cfmo  opanlng) 

(lat*  ( (naw- fiald-point-  (maka-naw-fiald-point-positiva-stap 

axposad-wall  fiald-point  stap-slza  ^2)) 
(cf-  (conflguratlon-f actor-for-ona-opanlng 

naw- fiald-point-  opanlng  axposad-wall) ) ) 

(cond  ( (>  cf-  (+  cfmo  .00001)) 

(saarch-for-max-mlnus 
axposad-wall  naw-flald-polnt- 
stap-slza  |32  cf-  opanlng)  ) 

(t  (altar-point  fiald-point  ;We  add  info  to  the  field  point. 

cf  cfmo ) 
fiald-point) ) ) ) 


5.  Determine  Area  to  Scan  for  Maximum 

We  now  wish  to  determine  the  area  to  scan  to  determine  the  maximum  configuration  factor 
in  some  region  of  an  exposed  wall.  We  note  that  the  radiation  pattern  due  to  a single  opening 
is  smooth  and  has  a single  peak.  The  radiation  from  more  than  one  opening  adds 
incoherently,  i.e.,  their  intensities  add.  Therefore,  a little  reflection  leads  one  to  conclude  that 
the  maximum  due  to  openings  will  lie  inside  the  area  determined  by  wrapping  a string  around 
pegs  placed  at  the  locations  of  the  individual  opeinings’  maximums.  However,  rather  than 
using  the  area  determined  with  the  pegs  and  string  we  will  determine  a rectangular  region 
that  lust  caputures  all  the  maximums.  In  some  cases,  such  as  for  openings,  there  will  be 
smaller  rectangular  regions  excuded  from  the  region  to  be  scannea. 
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We  begin  by  defining  a function  to  collect  the  location  of  all  the  maximums  of  the  individual 
openings.The  following  function  takes  as  input  the  liat-of ~«xposing-op«ning8  and  the 
«xpos«d-waii  and  returns  a list  of  the  maximum  point  for  the  openings.  The  slot  values  for 
these  points  includes  their  location  and 
confi^ration  factor. 


(defun  inaximum-of-an-open±ng-for-n-opening8 

(ll8t-of-e3Epo8±ng-°openi.ngs  expo8ed-waXl) 

(laap  ' li.8t  #'  (lambda  (opening) 

(Maxi miim-of-an-Openlng  opening  expoaed-wall) ) 
llst-of-expoalng-openlnga) ) 

We  will  have  use  for  the  function  aub-iiat  which  take  as  input  an  element  and  a 
ii8t-of  ”ii8t.  It  returns  a list  of  all  the  sub-lists  that  has  element  as  a member. 

(defun  8\ib-llat  (element  llat-of ~ll8t) 

(remove  element  llat-of-llst  :te8t-not  #' member)) 

We  need  two  functions.  The  first  determines  the  separations  between  the  projections  into  the 
xy  plane  of  a-point  and  each  point  in  a iiat-pointa. 

(defun  llst-xy-aeparatlons  (a-polnt  llat-polnta) 

(loop  for  point  In  llat-polnta 
collect 

(list  point  (xy-aeparatlon  a-polnt  point) ) ) ) 

The  second  function,  mln-aeparatlon-xy,  USeS  llat-xy-separatlons  tO  determine 
which  point  in  llat-polnta  which  has  the  smallest  xy-aeparatlon  from  a-polnt. 

(defun  mln-aeparatlon-xy  (a-polnt  llat-polnta) 

;; returns  the  point  nearest  a-polnt 
(let*  ( (llat-polnt-aep  (llat-xy-aeparatlona  a-polnt  llat-polnta) ) 

(llat-aep  (map  'Hat  #' aecond  llat-polnt-aep)) 

(mln-aep  (apply  'min  llat-aep))) 

(loop  for  polnt-aep  In  llat-polnt-aep 

when  (=  mln-aep  (aecond  polnt-aep) ) 
do 

(return  (first  polnt-aep) ) ) ) ) 


For  a region  which  has  no  openings  in  it  we  proceed  as  follows.  From  the  list  produced  by 
maximum-of-an-openlng-for-n-openlnga,  we  SelcCt  OUt  the  Smallest  and  the  largest 
values  for  z.  Then  we  use  xy-separation  to  determine  the  points  nearest  the  projection  of  the 
lower  left  and  upper  right  comers  of  the  exposed  region.  This  then  defines  a rectangular  area 
to  scan  that  will  include  the  maximum  of  the  openings.  We  used  the  rectangular  area  rather 
than  the  area  determined  by  wrapping  a string  around  the  pegs  because  it  was  easier  to 
determine. 

The  rectangular  regions  with  excluded  region  or  the  holey-rectangular  region  consists  of  a 
rectangular  region  with  rectangular  holes  cut  in  it.  If  a maximum  falls  into  an  opening,  we 
need  to  find  a replacement  maximum  for  this  errant  maximum.  This  replacement  maximum  will 
lie  on  the  bound^  of  the  the  opening  which  the  original  maximum  fell  into.  Thefore,  the 
region  we  scan  is  the  rectangular  region  as  determined  for  a nonholey  region  with  the 
scanning  area  expanded  to  include  the  opening. 

To  determine  if  a maximum  falls  in  an  opening,  we  need  a function  that  determines  whether 
a point  lies  in  the  interior  of  a rectangle.  The  function  point-in-r«ctangi4a  returns  the 
varticai-rectangi®  if  the  point  Hes  in  the  rectangle  and  nil  if  not. 


(defvm  point-in-ractangle  (vertical-rectangle  point) 

(let*  ( (11  (vertical-rectangle-lower-left  vertical-rectangle)  ) 
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(ur  (vax^±caJ.-rac^angla-upp«r- right  vartlcal-ractangla)  ) 

(xO  (point-x  11) ) 

(xl  (point-x  ur) ) 

(rO  (point- z 11) ) 

(zl  (polnt-z  ur) ) 

(yO  (polnt-y  11) ) 

(yl  (polnt-y  ur) ) 

(xu  (siax  xO  xl)) 

(xll  (min  xO  xl) ) 

(yu  (max  yO  yl) ) 

(yll  (min  yO  yl) ) 

(zu  (max  zO  zl) ) 

(zll  (min  zO  zl)) 

(xp  (polnt-x  point) ) 

(yp  (polnt-y  point) ) 

(zp  (polnt-z  point) ) ) 

(cond  ( (not  ( » xl  xO) ) 

(let*  ( (slope  (slope-of-vertlcal-rectangle  vertical-rectangle) ) 
(b  (-  yl  (*  slope  xl) ) ) ) 

;;  If  the  point  lies  on  the  boundary,  it  is  counted  as 
;;  outside  the  opening 

(cond  ( (and  (and  (>  zp  zll)  (<  zp  zu) ) 

(2uid  (>  xp  xll)  (<  xp  xu)  ) 

(=  (round-off  yp  4) 

(round-off  (+  (*  slope  xp)  b)  4) ) ) 
vertical - rectangle ) 

(t  nil)))) 

(t  (cond  ( (and  (and  (>  zp  zll)  (<  zp  zu) ) 

(auid  (>  yp  yll)  (<  yp  yu) ) ) vertical-rectangle) 
(t  nil)  ) ) ) ) ) 


Having  the  results  for  one  opening,  we  wish  to  determine  if  a point  lies  in  any  of  a list  of 
openings. 

(deftin  polnt-ln-an-openlng  (llst-of-openlngs  point) 

(loop  for  opening  In  llst-of-openlngs 

do  ;A  point  can  only  be  in  only  one  opening. 

(cond  ( (polnt-ln-rectangle  opening  point) 

(return  (polnt-ln-rectangle  opening  point) ) ) ) 
finally  (return  nil) ) ) 

Now  suppose  we  have  a list  of  points  and  a list  of  openings  and  we  want  a list  of  openings 
that  have  at  least  one  of  the  points  in  it. 

(defun  pointed-openings  (llst-of-openlngs  llst-of -points) 

(remove  nil  (loop  for  point  In  llst-of -points 

collect 

(polnt-ln-an-op>enlng  llst-of-openlngs  point)  ) ) ) 

If  we  have  a list  of  rectangles,  we  want  to  generate  a list  of  the  lower  left  and  upper  right 
comers. 

(defun  llst-of-comers  (llst-of-rectangles) 

(append 

(map  'list  #'  vertlcal-rectangle-lower-left  llst-of-rectangles) 

(map  'list  #'  vertical- rectangle-upper-right  llst-of-rectangles))) 

We  are  now  in  a position  to  determine  the  area  to  be  scanned  on  a wall  which  has  openings 
in  iL 


(deftin  area-to-acam-on-wall 
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(J.±afc~o£-«t3cposi.ng-opAn±ng8  «zposad-walJ.  H.8'b-o£-«zpo8ad-op>an±ngs) 
(l«t*  ( (l±8t-of-iftajE  ( nuixlmiim-of -an- opening- for— n-opanljiga 

l±8f-of-«zpo8lAg-op«ni^g8  azpoaad-wall) ) 
(li.8f-of-com«r8  (polntad-opanijiga  li.8t:-o£-«xpo8«d-opanlng8 

liat-of-max) ) 

;;If  iiat-of-com«ra  is  nil,  appand  docs  not  add  it  to  the  list. 
(li.st^-o£-max-plu8  (appand  li.af-o£-max  ld,sf-of-comara) ) 

(llaf-z  (map  ' li.8f  #'poinf-z  l±8f-of-max-plu8) ) 

(zl  (apply  'min  llaf-z) ) 

(z2  (apply  'max  lla^-z) ) 

(11  (wall-lowar-laf^  axpoaad-wall) ) 

(ur  (wall-uppar-rlghf  axpoaad-wall) ) 

(polnf-naar-11  (mln-aaparaflon-zy  11  llaf-of-max-plua) ) 
(poinf-naar-ur  (mln-aaparaflon-zy  \ir  llaf-of-max-plua)  ) 
(n-lowar-la£f  (maka -point:  :x  (polnf-x  point: -naar- 11) 

:y  (polnt-y  polnf-naar-11)  :z  zl) ) 
(n-uppar-right  (maka-poln't  :x  (polnt-x  polnt-naar-ur) 

:y  (polnt-y  polnt-naar-ur)  :z  z2))) 

(maka-vartlcal-ractangla 
: lowar-laft  n-lowar-laft 
:uppar-rlght  n-uppar- right) ) ) 


6.  Maximum  Configuration  Factor  for  each  subregion  of  the  exposed  wall. 

We  wish  to  determine  the  maximum  radiation  field  (or  configuration  factor)  for  each 
subregion  of  the  exposed  wall.  The  size  of  fields  we  are  interested  in  are  those  fields  that  can 
produce  piloted  ignition.  If  the  radiation  field  is  small  relative  to  that  required  for  piloted 
Ignition,  we  are  not  interested  in  it.  We  can  define  small  by  looking  for  the  most  vulnerable 
material  that  might  reasonably  be  exposed.  The  most  vulnerable  combustible  target  is  normally 
assumed  to  be  oven  dried  wood  [5].  If  the  maximum  possible  field  is  below  the  piloted 
ignition  level  for  this  material  then  it  is  assumed  that  no  exposure  problem  exists. 

To  take  advantage  of  this  assumption,  we  define  a global  constant 
*mini mum-pilot-ignition-lav«l*  tO  have  this  value  SO  we  Can 
use  it  as  the  default  value  for  the  piloted  ignition  level  for  materials. 


(dafvar  *minijmjm-pilot-ignition-l«v«l*  1.255 
"KlnimTim  pllohad  Ignlhlon  laval  for  ovan  drlad  wood  In  W./cm.  8q.") 

To  determine  the  maximum  configuration  factor  we  can  scan  the  entire  subregion  of  interest 
and  look  for  the  maximum.  However,  calculating  the  configuration  factor  for  all  these  points 
of  the  scan  takes  a lot  of  computer  computation  time.  Therefore,  it  is  desirable  to  hold  the 
number  of  points  we  scan  to  some  reasonable  minimum.  This  is  achieved  by  scanning  over 
the  area  in  steps  as  large  as  possible,  but  small  enough  so  we  don’t  make  a serious  error  in 
determining  the  maximum. 


(dafun  mln-8eap-8±za  (axpo8ad-wall  llsf-of -opanlngs) 
(loop  for  opanlng  In  li.8b-of-opani.ng8 
minimi ze 

(8bap-8iza  opening  axpoaad-wall) ) ) 


(dafun  nujnbar-8bapa  (m-8bap-aiza  ragion-bo-8can) 

(cond  ((or  (=  m-sbap-aiza  0)  ;Don't  want  to  divide  by  0. 

(<  (/  (varbical-racbangla-widbh  ragion-bo-acan) 

m-abap-aiza)  1))  1)  /Want  at  least  one  step. 

(b  (ro\ind  (/  (varbical-racbangla-widbh  ragion-bo-acan) 
m-abap-aiza) ) ) ) ) 


The  maximum  configuration  factor  due  to  a number  of  openings  is  always  less  than  the  sum 
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of  the  maximums  of  the  individual  openings.  We  will  have  use  for  this  sum  so  the  following 
function  computes  it 

(da£im  suin-o£-lnd±vi.du2Ll-c£  (l±8t-o£-op«iii.ngs  «3qpo8ad-waXl ) 

(apply  ' + (map  ' list:  #'  (lambda  (maz-polnb) 

(polnb-c£  max-poln't) ) 

(maatimum-of-an-opaning-for-n-opanlngs  ;list  of  nUlX  poifltS 

llsb-of-opaiilngs  «xpos«d-wall) ) ) } 

We  will  scan  the  wall  by  using  the  vertical  scan  to  handle  the  horizontal  scan.  If  a point 
falls  in  an  opening,  it  will  be  excluded. 

(defun  £lnd-nuuclmxim~con£lgurat:lon-£act:or-op«iilng 
(ll8b-o£-«3cpo8lng-opazilngs  azposad-wall 

&key  llsb-of-azposad-opanlngs 
(m±n±mum-pilot-ignition-lav«l“for-wall 
*in±nimim-pilot  - ignition- lavel  * ) 

;;the  default  value 

(sourco-flux  20) 
print-scan) 

(lat*  ( (tamp-maz-point  (maka-point  :x  0 :y  0 0 : cf  0)) 

(tamp-list-max-point  (list  tamp-max-point) ) ) 

(cond  ( (null  list-of-axposing-opanings) 

;;The  location  of  the  point  is  not  important  since  cf=0  which 
;;is  the  critical  information. 

tamp-list-max-point) 

((=1  (langth  list-o£-axposing-opanings) ) 

;;This  is  the  simple  case  of  one  opening.  Then  the  max  cf  is 
;;just  what  it  is  for  an  opening  & we  are  done. 

(list  (Maximum- o£-an-Opaning 

(£irst  li8t-o£-axposing-opaning8)  axposad-wall) ) ) 

;;Ifwe  have  more  than  one  opening  we  have  to  do  some  scanning  if 
; .there  is  a possibility  of  exceeding  threshold. 

(t  (lat  ((c£m  0) 

(s\jm-o£-c£s 

(stim-o£-individual-c£  list-o£-axposing-opanings 

axposad-wall) ) ) 

(cond  ((<  (*  1.25  sum-o£-c£s)  ;If  sum  of  all  moxs  is 

.’.’significantly  less  than  threshold,  we  don’t  need  to 
‘.’.calculate  cfbut  we  use  the  sum. 

(/  minimum-pilot-ignition-laval-£or-wall 
sourca-£lux) ) 

;;We  take  the  minimum-pilot-ignition-level-for-wall  and  divide  it 
;;by  source-flux,  which  has  a default  value  of  20  WIcm  sq.,  to  get 
.’.’minimum  configuration  factor. 

(print  'sum-o£“c£s) 

(satq  c£m  (round-o££  siam-o£-c£s  4)  ) 

;;The  following  is  not  the  location  of  the  max  or  its  value. 

; fit  doesn’t  matter  since  cf  is  below  threshold. 

(list  (altar-point  tamp-max-point  c£  c£m) ) ) 

;;  We  step  through  the  region  bounding  the  individual  maxs. 

The  vertical  scan  includes  the  horizontal  scan. 

(t  (scan-c£-vartically 

list-o£- axpo s ing - opanings  axposad-wall 
list-o£-axpoaad-opaninga  print-scan) ))))))) 


The  following  function  scans  the  horizontal  scanning  function  vertically. 


(dafun  scan-cf -vartically 

(list -of -expos ing- opening 8 exposed- wall 

list -of -exposed -openings  print -scan) 
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(l«t*  ( (r«g ion- to- scan 

(araa-to- scan-on- wall 

list-of-asqposing-opanings  «3cposad-wall 
list-o£-«xposad-op«nings) ) 

(zl  (point-z  (voirtical-ractangla-lowar-laft  ragion-to-acan)  ) ) 

(z2  (point-z  (vartical-ractangla-uppar-right  ragion-to-scan) ) ) 
(m-stap-siza  (xain-stap-siza  ascposad-wall 

list-o£-a3epoaing-opaninga)  ) 

(za  (point-z  (vartical-ractangla-lowa£’-la£t  ragion-to-scan)  ) ) 

(ys  (point-y  (vartical-ractangla-lowar-la£t  ragion-to-scem) ) ) 
(starting-point 

(maka-point  :z  zs  :y  ys  :z  zl) ) 

( initial - c£  ( con£iguration- £actor- £or-n- opanings 

starting-point  list-o£-azposing-opanings 
azposad-wall) ) 

(tamp-list-maz-point 
(list  (maka-point  :z  zs  :y  ys  :z  zl 

: c£  initial-c£) ) ) ) 

;;For  diagnostic  purposes  it  is  sometimes  of  interest  to  print  out 
;;the  area  to  be  scanned. 

(cond  (print -scan 

(£ormat  t "-%zl  = ■'2$,  yl  = ~2$,  zl  = ~2$"  zs  ys  zl) 

(format  t "~%zu  = ~2$,  yu  = -2$,  zu  = -2$" 

(point-z  (vartical-ractangla-uppar-right  ragion-to-sccui) ) 
(point-y  (vartical-ractangla-uppar-right  ragion-to-scan) ) 
z2)  ) ) 

;;Back  to  business. 

(loop  for  z from  zl  to  z2  by  m-stap-siza 
do 

(satq  tamp-list-maz-point 

(scan-cf -horizontally 

z list-of-azposing-opanings  azposad-wall 
list-of-azposad-opanings  tamp-list-maz-point 
m-stap-siza  ragion-to-scan  print-scan) ) 
finally  (raturn  tamp-list-maz-point) ) ) ) 

The  following  is  the  horizontal  scanning  function. 

(dafun  scan-cf-horizontaiiy  ;W e fix  z and  scan  over  xy . 

(z  list-of-azposing-opanings  azposad-wetll 

list-of-azposad-opanings  tamp-list-maz-point 
m-stap-siza  ragion-to-scan  print-scan) 

(lat*  ( (z  (point-z  (vartical-ractangla-lowar-laft  ragion-to-scan) ) ) 

(y  (point-y  (vartical-ractangla-lowar-laft  ragion-to-scan) ) ) 
(fiald-point  (maka-point  :z  z :y  y :z  z)) 

((32  (slopa-of-wall  azposad-wall)  ) 

(numbar-staps  (numbar-staps  m-stap-siza  ragion-to-scan) ) ) 

(loop  for  k from  0 to  numbar-staps 
do 

(altar-point  fiald-point 

z (+  z (*  k m-stap-siza  (cos  (32)  ) ) 
y (+  y Ic  m-stap-siza  (sin  (32)))) 

;;lf  field-point  is  in  an  opening,  we  skip  it. 

(cond  ( (point- in- an-opaning  list-of-azposad-opanings  fiald-point) ) 

;;Otherwise  we  process  it. 

(t  (lat 

( (naw-cf 

(round-off  (conf iguration- factor- for-n-opanings 

fiald-point  list-of-azposing-opanings 
azposad-wall)  4) ) 

(currant -max -cf  (point-cf  (first 

tamp-list-maz-point)  ) ) ) 
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;;If  we  make  print-scan  true,  we  get  a printout  of  the  scan. 

(cond  (prJLn'b-scan 
(foznoiat 
t 


"-%x  = ~2,0,6$,  y = ~2,0,6$,  r = ~2,0,6$,  n«w-cf 

a^ap-8i.za  = -’2,0,6$" 

(po±ne-z  flald-polnt:)  (point-y  fiald-polnt ) 

Id-point) 

naw-c£  m-stap-siza) ) ) 

(cond  ( (>  naw-cf  currant -max- cf)  ^replace  entire  Ust 
(aatq  tamp-liat-max-point 

(liat  (maka-point  :x  (point-x  fiald-point) 

:y  (point-y  flald-point) 

:z  (point-z  fiald-point) 

: cf  naw-cf) ) ) ) 

( (=  naw-cf  currant -max- cf)  ,'add  tO  Ust 
(aatq  tamp-liat-max-point 

(cona  (maka-point  :x  (point-x  fiald-point) 

•y  (point-y  fiald-point) 

:z  (point-z  fiald-point) 

:cf  naw-cf) 

tamp-liat-max-point) )))))) 
finally  (raturn  tamp-liat-max-point) ) ) ) 


-4,0,6$,  m- 
(point-z  fie 


To  convert  to  field  values  from  configuration  factors  we  multiply  configuration  factors  by 
the  source’s  flux.  This  is  such  a simple  thing  we  do  not  need  a function  to  compute  it. 


7.  Conclusions  and  Possible  Future  Work 

This  report  documents  a program  module  in  the  expert  system  EXPOSURE.  It  is  used  for 
calculating  the  maximum  electromagnetic  radiation  incident  upon  a building’s  wall  due  to  the 
burning  of  a neighboring  building.  It  is  assumed  that  the  burning  building  has  an  arbitrary 
number  of  rectangular  openings  emitting  radiation.  The  exposed  wall  can  be  considered  as 
being  composed  of  a number  of  rectangular  regions.  These  regions  may  be  openings  or 
regions  of  different  material.  The  results  of  calculating  the  configuration  factor  using  the 
function  configTiration-factor-for-on«-op«ning  with  analytical  expressions  for  three 
special  cases  showed  agreement  to  4 or  5 decimal  places  (see  Appendix  A).  This  is  just  little 
more  accuracy  that  the  expert  system  requires  so  that  is  desirable. 

This  module  also  allows  the  consideration  of  nonparallel  exposing  and  exposed  walls  and  the 
explicit  use  of  the  appropriate  wall  failure  mechanism.  In  the  case  of  comoustible  walls,  this 
means  using  the  piloted  ignition  threshold  for  the  combustible  material  rather  than  a very 
conservative  value  for  all  combustible  materials  such  as  is  done  by  NFPA  80A.  Also,  the 
effect  of  the  location  of  opening  is  taken  into  consideration.  Again  this  is  beyond  current 
practice.  These  features  improve  the  accuracy  of  the  analysis  of  the  safety  of  buildings  as  it 
relates  to  the  spread  of  fire  oetween  buildings  and  facilitates  innovative  design  and  use  of 
materials. 

Future  work  could  include  work  on  the  selection  of  the  step  size  used  in  scanning.  In  many 
real  world  problems  an  exposing  building  with  many  windows  exposes  a long  wall.  The 
present  program  could  prooably  be  made  faster  by  a more  intelligent  selection  of  the  initial 
step  size  and  increasing  the  step  size  in  regions  of  small  configuration  factors.  Also,  much 
more  effort  needs  to  be  spend  on  testing  the  program  before  it  could  be  considered  being 
used  in  a deployed  expert  system. 
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Appendix  A:  Configuration  Factor  Calculation  Compared  to  Results  of  Using  an  Analytic 
Expression 


In  this  appendix  we  will  compare  the  program’s  results  for  the  configuration  factor  of  a single 
opening  with  results  obtained  by  the  use  of  analytical  expressions  for  three  special  cases. 

Case  1. 

The  first  special  case  is  when  the  infinitesimal  element  is  parallel  to  the  opening  and  located 
on  a line  passing  through  a comer  of  the  opening  and  the  line  is  perpendicular  to  the  plane 
of  the  opening. 

The  analytical  expression  which  we  obtained  from  reference  8 page  823  #4,  is  expressed  in 
Lisp  code  in  the  function  cf-parall«l-corn«r. 

(d«£un  cf-parallttl-cornor  ( opanlng-wldth  opani.ng-halgh'k  aaparaeion) 

(loe*  ( (z  (/  opanJLng-wldeii  aaparaeion) ) 

(y  (/  op«iii.ng-*hai.ghe  saparaeion)  ) 

(sx  (sqrt  (+  1 (axpt  x 2)))) 

(ay  (sqrt  (+  1 («xpt  y 2)))) 

(tany  (aean  (/  y sx) ) ) 

(tanx  (a-ban  (/  x sy)  ) ) } 

(/  (+  (/  (*  tany)  sx)  (/  (*  y tanx)  sy)  ) (*  2 *single-p-pi*)  ) ) ) 

We  write  the  following  function  to  compare  the  program’s  result  with  that  obtained  by  using 

cf-parallol- comar. 

(dafun  program-c£-parallal-comar 

(opanij[\g -width  opanlng-halght  saparatlon) 

(lat*  ( (flald-polnt  (maka-polnt  :x  0 :y  saparatlon  :z  0)) 

;;  Field  point  Is  on  y axis. 

(llo  (maka-polnt  :x  opanlng-wldth  :y  0 ;z  0)) 

(uro  (maka-polnt  :x  0 :y  0 :z  opanlng-halght)) 

(opanlng  (maka-opanlng  : lowar-la£t  llo  :uppar-rlght  uro)) 

(llw  (maka-polnt  :x  -1  :y  saparatlon  :z  -1)) 

(urw  (maka-polnt  :x  1 :y  saparatlon  :z  1)) 

(axposad-wall  (maka-wall  :lowar-laft  llw  :uppar-rlght  urw))) 

(conf  Igruratlon- f actor- for-ona- opanlng 
flald-polnt  opanlng  axposad-wall) ) ) 

(dafun  slngla-comparlson-cf -parallal-comar 

(opanlng-wldth  opanlng-halght  saparatlon) 

(lat  ( (analytical 
(round-off 

(cf-parallal-comar  opanlng-wldth  opanlng-halght  saparatlon)  6)  ) 
(programs 

(program-cf -parallal-comar 

opanlng-wldth  opanlng-halght  separation) ) ) 

(format  t "~%~2$  ~6$  ~6$  ~6$  ~0$"  separation  analytical  programs 

(/  analytical  programs)  (*  (axpt  10  6)  (-  analytical  programs))))) 


(dafun  llst-comparlson-cf -parallal-comar  (llst-of-saparatlons ) 

(format  t "~%Y  analytical  program  analyt . /prog . dlffaranca  xl0^6") 
(loop  for  y In  llst-of-saparatlons 
do 

(slngla-comparlson-cf -parallal-comar  1 1 y)  ) ) 

If  we  mn  the  function 

(llst-comparlson-cf -parallal-comar  ' (.25  .5  1 2 3 4 5 6 7)) 
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with  the  indicated  arguments,  we  obtain: 


Y 

analytical 

program 

analy t . /prog . 

difference 

0.25 

0.237856 

0.237857 

0.999995 

-1. 

0.50 

0.207757 

0.207767 

0.999954 

-10. 

1.00 

0.138532 

0.138572 

0.999710 

-40. 

2.00 

0.059864 

0.059909 

0.999255 

-45. 

3.00 

0.030829 

0.030864 

0.998858 

-35. 

4.00 

0.018369 

0.018390 

0.998866 

-21. 

5.00 

0.012089 

0.012105 

0.998641 

-16. 

6.00 

0.008527 

0.008535 

0.999078 

-8. 

7.00 

0.006324 

0.006329 

0.999242 

-5. 

Case  2. 


Configuration  factor  perpendicular  over  one  comer 

The  second  special  case  is  when  the  infinitesimal  element  is  perpendicular  to  the  opening  and 
located  on  a fine  passing  through  a comer  of  the  opening  and  the  fine  is  perpendicular  to  the 
plane  of  the  opening. 

The  analytical  expression  which  we  obtained  from  reference  8 page  823  #6,  is  expressed  in 
Lisp  code  in  the  function  cf -p«rp«ndi.cular-eomar . 

(daf\m  c£-p«rp«nd±cular~cornar 

(op«ni.ng-wi.dt;h  opani.ng-halght:  separation} 

(let*  ( (x  (/  opening- width  opening-height) ) 

(y  (/  separation  opening-height) ) 

(sxy  (sqrt  (+  (expt  x 2)  (expt  y 2)))) 

(tany  (atan  (/  1 y) ) ) 

(tanxy  (atan  (/  1 sxy) ) ) ) 

(/  (-  tany  (/  (*  y tanxy)  sxy) ) (*  2 *slngle-p-pl*) ) ) ) 


(defun  program-cf-perpendlcular-comer 

(opening-width  opening-height  separation) 

(let*  ((field-point  (make-point  ;x  0 :y  separation  :z  0)) 

Field  point  is  on  y axis. 

(llo  (make-point  :x  opening-width  :y  0 :z  0)) 

(uro  (make-point  ;x  0 :y  0 :z  opening-height)) 

(opening  (make-opening  ; lower-left  llo  : upper-right  uro)) 
(llw  (make-point  :x  0 :y  separation  :z  -1)) 

(urw  (make-point  :x  0 :y  (-i’  separation  1)  :z  1)) 
(exposed-wall  (make-wall  : lower-left  llw  : upper-right  urw))) 
(configuration- factor- for-one-openlng 
field-point  opening  exposed- wall) ) ) 


(defun  single-comparison-cf-perpendicular-comer 
(opening- width  opening-height  separation) 

(let  ( (analytical 
(round-off 

( cf  -perpend!  cular  - comer  opening  - width 

opening-height  separation)  6) ) 

(programs 

(program-cf-perpendicular-corner 
opening -width  opening- height  separation) ) ) 

(format  t "‘>'%~2$  -6$  ~6$  ~6$  ~0$"  separation  analytical  programs 

(/  analytical  programs)  (*  (expt  10  6)  (-  analytical  programs))))) 


(defun  list- comparison -cf -perpendicular- comer 
(list-of -separations ) 


S : >smith> radiation . lisp. 33 


9/24/90  13:50:41  Page  23 


(format  t '•~%Y  analytical  program  analyt . /prog . dlfforanco  xlO'^S") 
(loop  for  y in  liat-of-aoparationa 
do 

(aingla-cou^ariaon-cf-porpondicular-comar  1 1 y) ) ) 

If  we  run  the  function 

(liat-compariaon-cf-porpandicular-comar  ' (,5  1 2 3)) 

with  the  indicated  arguments,  we  obtain: 


Y 

analytical 

program 

analyt . /prog. 

difference 

0.50 

0.124269 

0.124326 

0.999541 

-57. 

1.00 

0.055734 

0.055789 

0.999010 

-55. 

2.00 

0.013928 

0.013952 

0.998309 

-24. 

3.00 

0.004964 

0.004976 

0.997573 

-12. 

Case  3. 

Configuration  factor  angle  theta  over  one  comer 

The  third  special  case  is  when  the  infinitesimal  element  is  in  a plane  that  has  an  angle  of 
intersection  with  the  plane  of  the  opening  of  theta.  The  separation  is  the  distance  from  the 
comer  of  the  opening  where  the  two  planes  intersect  to  the  infinitesimal  element. 

Reference  "Fire-Safe  Stmctural  Steel  A Design  Guide"  American  Iron  and  Steel  Institute, 
(1979)  page  24 


(defim  cf-bhata-comer  (opaning-width  opaning-haighb  saparabion  bhaba) 
(lab*  ( (x  (/  opaning-widbh  aaparabion) ) 

(y  (/  opaning-haighb  saparabion) ) 

(yc  (*  y (coa  bhaba) ) ) 

(ay  (aqrb  (-  (+  1 (axpb  y 2) ) (*  2 yc)  ) ) ) 

(xa  (aqrb  (+  (axpb  x 2)  (axpb  (ain  bhaba)  2) ) ) ) ) 

(/  (+  (aban  x)  (/  (*  (-  yc  1)  (aban  (/  x ay)  ) ) ay) 

(*  (*  X (coa  bhaba) ) (/  (+  (aban  (/  (-  y (coa  bhaba) ) xa) ) 

(aban  (/  (coa  bhaba)  xa) ) ) xa) ) ) 

(*  2 *aingla~p-pi*) ) ) ) 


(dafun  program-cf-bhaba-cornar 

(op>anlng-widbh  opaning-haighb  aaparabion  bhaba) 

(lab*  ( (fiald-poinb  (maka-poinb  :x  (*  aaparabion  (coa  bhaba)) 

:y  (*  aaparabion  (ain  bhaba)) 

:z  0)) 

(llo  (maka-poinb  :x  opaning-widbh  :y  0 :z  0)) 

(uro  (maka-poinb  :x  0 :y  0 :z  opaning-haighb)) 

(opaning  (maka-opaning  :lowar-lafb  llo  :uppar-righb  uro)) 
(llw  (maka-poinb  :x  (*  aaparabion  (coa  bhaba)) 

:y  (*  aaparabion  (ain  bhaba)) 

-D) 

(urw  (maka-poinb  :x  (*  (+  aaparabion  1)  (coa  bhaba)) 

:y  (*  (+  aaparabion  1)  (ain  bhaba)) 

:r  1)) 

(axpoaad-wall  (maka-wall  :lowar-lafb  llw  :uppar-righb  urw))) 
(conf igurabion-f acbor- for- ona- opaning 
fiald-poinb  opaning  axpoaad-wall) ) ) 


(dafun  aingla-compariaon-cf -bhaba-cornar 

(opaning-widbh  opaning-haighb  aaparabion  bhaba) 
(lab  ( (analybical 
(round-off 

(cf -bhaba-comar  opaning-widbh 
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op«ni.ng-halgh1:  separation 

theta)  6)  ) 

(programs 

(program-cf -theta-comer 

opening -width  opening-height  separation  theta) ) ) 

(format  t "-%~2$  ~3$  -6$  ~6$  ~6$  ~0$" 

separation  theta  analytical  programs 

(/  analytical  programs)  (*  (ezpt  10  6) (-  analytical  programs) ) ) ) ) 

(defnn  llst-comparlson-cf-theta-comer 

(llst-of-separatlons  llst-of -thetas ) 

(format  t "•'%Y  theta  analytical  program  analyt . /prog . difference  3e10‘^6”) 

(loop  for  xy  In  llst-of-separatlons 
do 

(loop  for  theta  In  llst-of-thetas 
do 

(slngle-comparlson-cf-theta-comer  1 1 xy  theta) ) ) ) 

If  we  run  the  function  (llst-con^arlson-cf-theta-comer  '(12  3)  (list  (/ 

*8lngle-p-pl*  2)  3 (/  *slngle-p-pl*  4))  (*  3.5  (/  *slngle-p-pl*  4) ) (* 

3.9  (/  *singie-p-pi*  4) ) ) ) with  the  indicated  arguments,  we  obtain: 


Y 

theta 

analytical 

program 

analyt , /prog . 

difference 

1.00 

1.571 

0.055734 

0.055789 

0.999010 

-55. 

1.00 

2.356 

0.013020 

0.013042 

0.998303 

-22. 

1.00 

2.749 

0.003201 

0.003207 

0.997986 

-6. 

1.00 

3.063 

0.000127 

0.000128 

0 . 995048 

-1 . 

2.00 

1.571 

0.013928 

0.013952 

0.998309 

-24. 

2.00 

2.356 

0.003797 

0.003810 

0 . 996505 

-13. 

2.00 

2.749 

0.000961 

0.000965 

0 . 996236 

-4  . 

2.00 

3.063 

0.000039 

0.000039 

1.007274 

0. 

3.00 

1.571 

0.004964 

0.004976 

0.997573 

-12. 

3.00 

2.356 

0.001542 

0.001552 

0.993681 

-10. 

3.00 

2 . 749 

0.000401 

0.000403 

0.993925 

-2  . 

3.00 

3.063 

0.000016 

0.000016 

0 . 980120 

-0  . 
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Appendix  B:  Results  for  Some  Simple  Buildings 

We  present  here  the  results  of  applying  the  results  of  this  file  to  three  cases  that  will  test 
the  performance  of  this  program.  It  is  not  a comprehensive  test,  but  one  that  should  exercise 
the  program’s  main  features. 

We  define  a function  that  will  take  as  its  first  argument  a list  of  the  list  of  lower  left  and 
upper  right  comers  of  the  openings.  The  next  six  arguments  are  the  coordinates  of  the  lower 
left  and  upper  right  comers  of  the  exposed  area, 

(daf un  'bast-  fi.nd-mazlmum-conflgura'blon-factor-opani.ng 
(list-of-paJLrs-of-cornars-of-opanings 

x3  y3  z3  x4  y4  z4  ;The  Corner  Coordinates  of  the  exposed  area. 

4kay  print-scan)  ;Scan  pattern  produced  upon  request. 

(lat* 

( (list-of-azposing-opanings 
(map  ' list 

#'  (lambda  (list-o£-pair-o£-comars) 

(lat  ((11  (maka-point  :z  (first  list-of-pair-o£-comars) 

:y  (sacond  list-o£-pair-o£-eomars) 

: z (third  list-o£-pair-o£-comars)  ) ) 

(ur  (maka-point  :x  (fourth  list-of-pair-of-comars) 

:y  (fifth  list-of-pair-of-comars) 

:z  (sixth  list-of-pair-of-comars)))) 
(maka-opaning  :lowar-la£t  11  :uppar-right  ur) ) ) 
lis t-of -pairs -of -comar s -of -opanings)  ) 

(llw  (maka-point  :x  x3 
:y  y3 
:z  z3)) 

(urw  (maka-point  :x  x4 
:y  y4 
: z z4) ) 

(axposad-wall  (maka-wall  :lowar-la£t  llw  :upp>ar-right  urw)  ) ) 

( find-maximum-configuration- factor- opaning 

list-of-axposing-opanings  axposad-wed.1  : print-scan  print-scan) ) ) 

We  will  test  three  cases.  In  all  cases  we  have  two  exposing  windows  either  in  the  xz  or  yz 
plane.  While  the  windows  are  the  same  size,  they  are  set  m at  different  heights.  For  the 
first  two  cases  the  exposed  wall  is  parallel  to  the  exposing  wall.  In  the  last  case,  it  is  at  a 
small  angle  to  the  exposing  wall. 

Our  first  test  examines  the  time  need  to  compute  the  results. 

(da£\m  tast-casas  () 

(values  (time  (test- find-maxi mum-configuration- factor-opening 

'((011077)  (0  9 7 0 15  13))  20  0 0 20  15  15)) 

(terpri) 

;; Interchanging  x and  y,  we  get 

(time  (teat- find-maxi mum-configuration- factor-opening 

'((101707)  (9  0 7 15  0 13))  0 20  0 15  20  15) ) 

(tarpri) 

;; Setting  the  exposed  wall  at  a small  angle  to  the  exposing  wall. 

(time  (teat- find-maximum-configuration- factor-opening 

'((101707)  (9  0 7 15  0 13))  0 18  0 15  22  15) )) ) 

The  results  of  running  test-cases  follows: 

(test-cases ) 

Evaluation  of  (TEST-FIND-MAXIMUM-CONFIGURATION-FACTOR-OPENING  ' # 20  0 0 . . . ) 
took  8.441768  seconds  of  elapsed  time  including: 

0.112  seconds  processing  sequence  breaks, 

0.360  seconds  in  the  storage  system  (including  0.052  seconds  waiting  for  pages): 
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0.179  seconds  processing  313  page  faults  including  4 fetches, 

0.171  seconds  in  creating  and  destroying  pages,  and 
0.010  seconds  in  miscellaneous  storage  system  tasks. 

The  garbage  collector  has  flipped;  so  no  consing  was  measured. 

Evaluation  of  (TEST-FIND-MAXIMUM-CONFIGURATION-FACTOR-OPENING  ' # 0 20  0 . . . ) 
took  6.812716  seconds  of  elapsed  time  including: 

0.137  seconds  processing  sequence  breaks, 

0.531  seconds  in  the  storage  system  (including  0.135  seconds  waiting  for  pages): 
0.312  seconds  processing  330  page  faults  including  21  fetches, 

0.159  seconds  in  creating  and  destroying  pages,  and 
0.060  seconds  in  miscellaneous  storage  system  tasks. 

235  list,  40,365  structure,  12,325  stack  words  consed  in  WORKING-STORAGE-AREA. 
Evaluation  of  (TEST-FIND-MAXIMUM-CONFIGURATION-FACTOR-OPENING  '#  0 18  0 ...) 
took  8.779506  seconds  of  elapsed  time  including: 

0.173  seconds  processing  sequence  breaks, 

0.587  seconds  in  the  storage  system  (including  0.082  seconds  waiting  for  pages); 
0.220  seconds  processing  391  page  faults  including  6 fetches, 

0.231  seconds  in  creating  and  destroying  pages,  and 
0.135  seconds  in  miscellaneous  storage  system  tasks. 

221  list,  44,207  structure,  18,121  stack  words  consed  in  WORKING- STORAGE -AREA. 
(#S(POINT  :X  20.0 

:Y  6.828427 
:Z  6.828427 
:CF  247/5000) ) 


NIL 

(#S(POINT  :X  6.828427 
:Y  20.0 
:Z  6.828427 
:CF  247/5000) ) 


NIL 

(#S(POINT  :X  2.5196323 
:Y  18.671902 
:Z  6.607681 
:CF  513/10000) ) 


The  second  shows  the  area  to  be  scanned  and  the  actual  scanned  points. 

(defun  test-cases-2  () 

(values  (tesb-f  ind-maxi, mum- configuration -fact  or-openlng 

'((011077)  (0  9 7 0 15  13) ) 20  0 0 20  15  15  :prlnt-scan  t) 

(terprl) 

// Interchanging  x and  y,  we  get 

(test-flnd-marimum-conflguratlon-f actor-opening 

'((1  0 1 7 0 7)  (9  0 7 15  0 13)  ) 0 20  0 15  20  15  :prlnt-scan  t) 

(terprl) 

//Setting  the  exposed  wall  at  a small  angle  to  the  exposing  wall. 

(test-flnd-maximum-conflguratlon- factor-opening 
'((101707)  (9  0 7 15  0 13)) 

0 18  0 15  22  15  :prlnt-8can  t) ) ) 

The  results  of  running  test-cases-2  follows: 

(test-cases-2 ) 


xl  = 

20 . 00, 

yl  = 

O 

O 

zl 

= 4 

.00 

XU  = 

20 . 00, 

yu  = 

12.00, 

zu  = 

10.00 

X = 

20.00, 

y = 

4.00, 

z 

= 

4.00, 

new-cf 

.0460, 

m- step- size 

2 . 83 

X = 

20.00, 

y = 

6.83, 

z 

= 

4.00, 

new-cf 

. 0479, 

m-step-slze 

= 

2.83 

X = 

20.00, 

y = 

9.66, 

z 

= 

4.00, 

new-cf 

= 

. 0471, 

m-step-slze 

= 

2.83 

X = 

20.00, 

y = 

12.49, 

z 

= 

4 . 00, 

new-cf 

SE 

. 0439, 

m-step-slze 

= 

2 . 83 

X = 

20.00, 

y = 

4.00, 

z 

= 

6.83, 

new-cf 

= 

. 0470, 

m-step-slze 

= 

2.83 

X = 

20.00, 

y = 

6.83, 

z 

= 

6.83, 

new-cf 

. 0494, 

m-step-slze 

= 

2 .83 

X = 

20 . 00, 

y = 

9.66, 

z 

= 

6. 83, 

new-cf 

= 

. 0491, 

m-step-slze 

= 

2 . 83 
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z 

= 

20.00, 

y = 

12.49, 

z = 

6.83, 

naw-c£ 

3 

. 0462, 

ni-s£ap-slza 

3 

2.83 

z 

= 

20.00, 

y = 

4.00, 

z = 

9.66, 

n«w-cf 

= 

. 0450, 

m-a'tap-slza 

3 

2 . 83 

z 

=s 

20.00, 

y = 

6,83, 

z = 

9,66, 

n«w-cf 

= 

.0478, 

m-s£ep-alza 

= 

2.83 

z 

= 

20.00, 

y = 

9.66, 

z = 

9.66, 

n«w-cf 

= 

.0481, 

m-abap-alza 

3 

2.83 

z 

20.00, 

y = 

12.49, 

z = 

9.66, 

n«w-c£ 

3 

. 0457, 

m-at:ap-alza 

3 

2.83 

zJ 



4.00, 

yl  = 

20.00, 

zl  = 

4 

.00 

zu  = 

12.00, 

yu  = 

= 20.00, 

zu  = 

10.00 

z 

s: 

6.83, 

y = 

20.00, 

z = 

4.00, 

naw-c£ 

3 

. 0479, 

m-abap-alza 

= 

2.83 

z 

= 

9.66, 

y = 

20.00, 

z = 

4.00, 

nttw-c£ 

=* 

.0471, 

m-  a'bap-  al  za 

3 

2.83 

z 

12.49, 

y = 

20.00, 

z = 

4.00, 

UttW-Cf 

3 

. 0439, 

m-atap-alza 

3 

2 . 83 

z 

= 

4.00, 

y = 

20.00, 

z = 

6.83, 

now-c£ 

3 

.0470, 

m-a£ap-a±za 

3 

2.83 

z 

= 

6.83, 

y = 

20,00, 

z = 

6.83, 

naw-c£ 

= 

.0494, 

m-a£ap-alza 

3 

2.83 

z 

=s 

9.66, 

y = 

20.00, 

z = 

6.83, 

naw-c£ 

3 

. 0491, 

in-a£ap-alza 

3 

2.83 

z 

3 

12.49, 

y = 

20.00, 

z = 

6.83, 

n«w-c£ 

3 

.0462, 

m-a£ap~alza 

3 

2.83 

z 

= 

4.00, 

y = 

20.00, 

z =* 

9.66, 

n«w-c£ 

3 

.0450, 

m^-atap-alza 

3 

2.83 

z 

3 

6.83, 

y = 

20.00, 

Z 

9.66, 

n«w-c£ 

3 

.0478, 

m-a£ap-alza 

3 

2.83 

z 

= 

9.66, 

y * 

20.00, 

z = 

9.66, 

naw-c£ 

3 

. 0481, 

m- a'bap- alza 

3 

2.83 

z 

= 

12.49, 

y = 

20.00, 

z = 

9.66, 

n«w-c£ 

3 

.0457, 

m-abap-alza 

3 

2.83 

zl 

0.00, 

yl  = 

18.00, 

zl  = 

4 

.00 

3£U  = 

6.72, 

yu  = 

19.79, 

zu  = 

10.00 

z 

= 

.00, 

y = 

18.00, 

z = 

4.00, 

naw-c£ 

3 

. 0485, 

m-abap-alza 

3 

2 . 61 

z 

= 

2.52, 

y = 

18.67, 

z = 

4.00, 

n«w-c£ 

3 

. 0503, 

m-abap-alza 

3 

2.61 

z 

= 

5.04, 

y = 

19.34, 

z = 

4.00, 

n«w-c£ 

3 

. 0495, 

m-abap-alza 

3 

2.61 

z 

= 

7.56, 

y = 

20.02, 

z = 

4 . 00, 

n«w-c£ 

3 

. 0463, 

m-abap-alza 

3 

2 . 61 

z 

s= 

. 00, 

y = 

18.00, 

z = 

6.61, 

now-c£ 

3 

. 0490, 

m-abap-alza 

3 

2.61 

z 

= 

2.52, 

y = 

18.67, 

z = 

6.61, 

naw-c£ 

3 

.0513, 

m-abap-alza 

3 

2.61 

z 

= 

5.04, 

y = 

19.34, 

z = 

6.61, 

n«w-c£ 

3 

. 0509, 

m-abap-alza 

3 

2.61 

z 

= 

7.56, 

y = 

20.02, 

z = 

6.61, 

naw-c£ 

3 

. 0481, 

m-abap-alza 

3 

2 . 61 

z 

3 

.00, 

y = 

18.00, 

z — 

9.22, 

naw-c£ 

3 

. 0467, 

m-abap-alza 

3 

2.61 

z 

= 

2.52, 

y = 

18.67, 

z = 

9.22, 

n«w-c£ 

3 

.0494, 

m-abap-alza 

3 

2 . 61 

z 

3 

5.04, 

y = 

19.34, 

z = 

9.22, 

n«w-c£ 

3 

.0495, 

m-abap-alza 

3 

2 . 61 

z 

= 

7.56, 

y = 

20.02, 

z = 

9.22, 

n«w-c£ 

3 

. 0473, 

m-abap-alza 

3 

2 . 61 

(#S  (POINT  ; 

X 20. 

0 

:Y  6.828427 
:Z  6.828427 
:CF  247/5000) ) 

NIL 

(#S(POINT  :X  6.828427 
:Y  20.0 
:Z  6.828427 
;CF  247/5000) ) 

NIL 

(#S(POINT  :X  2.5196323 
:Y  18.671902 
:Z  6.607681 
:CF  513/10000) ) 
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